arXiv:1503.08330vl [math-ph] 28 Mar 2015 


Resolution of Chern-Simons-Higgs Vortex Equations 


Xiaosen Han 

Institute of Contemporary Mathematics 
School of Mathematics 
Henan University 
Kaifeng, Henan 475004, PR China 

Chang-Shou Lin 

Taida Institute for Mathematical Sciences 
National Taiwan University 
Taipei, Taiwan 10617, ROC 

Yisong Yang 

Department of Mathematics 
Polytechnic School of Engineering 
New York University 
Brooklyn, New York 11201, USA 
NYU-ECNU Institute of Mathematical Sciences 
New York University - Shanghai 
3663 North Zhongshan Road, Shanghai 200062, PR China 


Abstract 

It is well known that the presence of multiple constraints of non-Abelian relativisitic 
Chern-Simons-Higgs vortex equations makes it difficult to develop an existence theory 
when the underlying Cartan matrix K of the equations is that of a general simple Lie 
algebra and the strongest result in the literature so far is when the Cartan subalgebra 
is of dimension 2. In this paper we overcome this difficulty by implicitly resolving 
the multiple constraints using a degree-theorem argument, utilizing a key positivity 
property of the inverse of the Cartan matrix deduced in an earlier work of Lusztig 
and Tits, which enables a process that converts the equality constraints to inequality 
constraints in the variational formalism. Thus this work establishes a general existence 
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theorem which settles a long-standing open problem in the field regarding the general 
solvability of the equations. 

Keywords: Chern-Simons-Higgs vortex equations, non-Abelian gauge theory, Lie 
algebras, Cartan matrices, constraints, calculus of variations, minimization. 
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1 Introdcution 

The main result of this paper is a proof of a general existence theorem for the doubly-periodic 
solutions of the relativistic non-Abelian Chern-Simons-Higgs vortex equations whose Cartan 
matrix K is that of an arbitrary simple Lie algebra. The motivation of our work originates 
from theoretical physics and development of new methods of mathematical analysis to tackle 
systems of nonlinear equations of physical interests. Below we begin with a description of 
some of the field-theory aspects of our study. 

Field-theoretical origins. In 1957, Abrikosov [I] predicted in the context of the Ginzburg- 
Landau theory [18113011551163] that vortices of a planar lattice structure may appear in super¬ 
conductors of the second type. Due to the complexity of the problem, a complete mathemat¬ 
ical understanding of such Abrikosov vortices has not been obtained yet, although in an over 
simplified situation, known as the Bogomornyi-Prasad-Sommerficld (BPS) limit [9j[53], an 
existence and uniqueness theorem is proved under a necessary and sufficient condition m 
where the Abrikosov lattice structure is realized through gauge-periodic boundary conditions 
conceptualized and formulated by ’t Hooft [63] so that the governing elliptic equation is 
defined over a doubly periodic planar domain. Since early 1980, there has been a growing in¬ 
terest in accommodating electrically charged vortices in condensed-matter physics (2714291169] 
such that the presence of the Chern-Simons dynamics becomes imperative due to the Julia- 
Zee theorem which says that vortices in the classical Yang-Mills-Higgs models can only 
carry magnetism. For a recent review on these subjects and related literature, see (72j- It is 
relevant to describe such a non-go theorem in the context of the simplest Yang-Mills-Higgs 
model here. 

The Julia-Zee theorem. Consider the classical Yang-Mills-Higgs (YMH) model in the 

adjoint representation of SU (2) over the Minkowski spacetime M 2,1 of signature (-1-) and 

use /a, v — 0,1,2 to denote the spacetime coordinate indices. Let A p and (j) be the gauge 
and Higgs fields written as isovectors, respectively. Then the action density is 

-2" = - jF- • F„„ + 1 D“4 ■ D,4> - - l) 2 , (1-1) 

where the field strength tensor F /1L/ is defined by F pu = d p A u — d v A p + eA p x A u and the 
gauge-covariant derivative D p is defined by D p cj) = d p (j) + eA p x qb. Here A, e are positive 
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coupling constants. The Euler-Lagrange equations of (II.ip in the static limit are 


( 1 . 2 ) 

(1.3) 

(1.4) 


AA0 = —e 2 (| A o | 2 0 — [4> ■ A 0 ]A 0 ) + A0(|0| 2 — 1), 

DjFij = e(A 0 x Fj 0 + f x Drf), 

AA 0 = —e((9j[Aj x A 0 ] + A* x [<9jA 0 + eAx A 0 ]) + e 2 f x (A 0 x 0), 

over the space R 2 , among which (11.41) is also the Gauss law constraint. The two-dimensional 

nature of these equations implies that their solutions may be interpreted as the YMH vortices. 
On the other hand, the associated energy density, or Hamiltonian, of (II. ip may be calculated 
to be 

Jf = -|Fjj| 2 + — |(9 ,;Aq + eA i x A 0 | 2 + A0| 2 + -e 2 |A 0 x f | 2 + — (\f>\ 2 — l) 2 , (1.5) 

which is positive definite. In [59] it is proved that under the finite-energy condition E = 
f R 2 J^dx < oo, any solution of (11.21) (1 1.4 f) must have a trivial temporal component for its 
gauge field: 

A 0 = 0. (1.6) 

This result, first observed in [4U] for the Abelian Higgs model, later referred to as the Julia- 
Zee theorem [59], has an important physical implication. Recall that, based on consideration 
on the manner of interactions, ’t Hooft [62] proposed that the electromagnetic field F )W in 
the YMH theory is to be defined by the formula 

F P* = - D u (p). (1.7) 

m e|0| 3 

Hence, inserting (11.61) into (11.71) . we obtain F 0 i = = 0, which indicates that there is 

no induced electric field in the model. In other words, the YHM vortices, Abelian or non- 
Abelian, are purely magnetic and electrically neutral. 

Emergence of the Chern-Simons type models. The development of theoretical physics 
requires the presence of both electrically and magnetically charged vortices, also called dyons, 
since such dually charged vortices have found applications in a wide range of areas including 
high-temperature superconductivity [43j,[4?], optics [8], the Bose-Einstein condensates m 
|42| . the quantum Hall effect [58], and superfluids HZHSHEZ]. Thanks to the studies of 
Jackiw-Templeton [38], Schonfeld T6| . Deser-Jackiw-Templeton 17*, 201- Paul-Khare [ 52] . 
de Vega-Schaposnik [65y66], and Kumar-Khare ra. starting from the early 1980’s, it has 
become accepted that, in order to accommodate electrically charged vortices, one needs to 
introduce into the action Lagrangian a Chern-Simons topological term [15, 16]. which has 
also become a central structure in anyon physics [29[[68j[69]. On the other hand, despite of 
the importance of electrically charged vortices with the added Chern-Simons dynamics, it 
has been a difficult issue until rather recently [T2j to construct finite-energy solutions of the 
field equations because of the indefiniteness of the action functional as a consequence of the 
Minkowski signature of spacetime and the presence of electricity. In 1990, it came as a fortune 
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that Hong, Kim, and Pac [35] and Jackiw and Weinberg |39| showed that when one uses 
only the Chern-Simons term and switches off the usual Maxwell term in the Abelian Higgs 
model one can achieve a BPS structure and thus arrive at a dramatic simplification of the 
governing equations. Subsequently, the ideas of [35j,:39] were extended to non-Abelian gauge 
held theory models and a wealth of highly interesting systems of nonlinear elliptic equations of 
rich structures governing non-Abelian Chern-Simons-Higgs vortices was unearthed [22H21]- 
More recently, these ideas have also been further developed in supersymmetric gauge held 
theory in the context of the Bagger-Lambert-Gustavsson (BLG) model j4H7l[Tl'j25,.l3T] and 
the Aharony-Bergman-Jafferis-Maldacena (ABJM) model [2111411331144] which have been the 
focus of numerous activities in contemporary held-theoretical physics. 

The present work is a complete resolution of the most general relativistic Chern-Simons- 
Higgs vortex equations defined over a doubly periodic planar domain with the Cartan matrix 
of an arbitrary simple Lie algebra. In the next section, we describe the vortex equations and 
some of the key technical issues, including methodology. In the section that follows, we state 
our main existence theorem. In the subsequent sections, we prove this theorem. 

2 Vortex equations, technical issues, and methodology 

Let K = ( Kij ) be the Cartan matrix of a finite-dimensional semisimple Lie algebra L. We 
are interested in the relativistic Chern-Simons-Higgs vortex equations [22h24LfTOllTI] of the 
form 

( n n n \ Ni 

y y K kj K ji e Uj e Uk - y Kjie Uj J + 47t y S Pij (x), i = 1,..., n, (2.1) 

j =1 k =1 j =1 / j=l 

where n > 1 is the rank of L which is the dimension of the Cartan subalgebra of L , 5 p 
denotes the Dirac measure concentrated at the point p, X > 0 is a coupling constant, and the 
equations are considered over a doubly periodic domain f! resembling a lattice cell housing 
a distribution of point vortices located at p t j , j = 1,..., N tl i = 1,..., n. For an existence 
theory the ultimate goal is to obtain conditions under which (12.ip allows or fails to allow 
a solution. In order to see the technical difficulties of the problem, we take the beginning 
situation n — 1 as an illustration for which the underlying gauge group may be either 17(1) 
or SU (2) which is of fundamental importance in applications, such that (12. ip takes the scalar 
form 

N 

Au = Xe u (e u — 1) + An y^ 5 Pj (x), x e fl (2.2) 

3=1 

Let «o be doubly periodic modulo Q and satisfy An 0 = — + 4x Vy=i VV)- Then 

u = uq + v recasts (12.2jl into 

Ajr N 

Av = \e U0+v (e U0+v - 1) + ——, (2.3) 

|iz| 
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which is the Euler-Lagrange equation of the action functional 


S(v) 




- Ae“ 0+W + 


47 tN 1 

¥j 


dx. 


(2.4) 


It is easily seeing by taking constants as test functions that (12.4p is not bounded from below 
over the space of doubly periodic functions. In order to tackle such a difficulty, it is to take 
into account of a natural constraint arising from integrating 12.31) . That is, 



|fi| 4vriV 
1 ~ 


(2.5) 


which indicates that there fails to permit a solution when A < 1 ^^ v . At a first glance, it 
may be attempting to believe that a solution may be obtained by minimizing (12.4p subject 
to (12.51) . Unfortunately, there arises a Lagrange multiplier issue which prohibits a minimiza¬ 
tion process under the equality constraint (12.5p . To overcome this issue, we may take the 
decomposition v = c + w where c6t and f Q w dx = 0 and rewrite the constraint (12.5p as 

r r A'jr ]\r 

e 2c / e 2u 0+ 2w dx _ Q c e u 0 +w ^ + = Q) ( 2 , 6 ) 

Jn Jo. A 

which becomes a solvable quadratic equation in £ = e c if and only if the discriminant of (12.6ft 
stays nonnegative, 

(^J e uo+w dxj - J e 2uo+2w dx > 0, (2.7) 

so that c may be represented as (say) 


c = In 


3 lt0+u> 



e «o+™ dx — 


lQnN 

A 


e 2uo+2w dx -hi 2 e 2no+2«, da; 


( 2 . 8 ) 

Then it may be shown that a solution of (12.3p can be obtained by minimizing the action 
functional 

I{w) = J |^|VH 2 + ^e 2uo+2c+2w - Ae“ 0+C+U, | dx + AttNc, (2.9) 

descending from (12.41) . subject to the inequality constraint (12.71) . with c given by (12. 8p . when 
A is sufficiently large so that a minimizer occurs in the exterior of the constraint class which 
rules out the Lagrange multiplier problem mentioned earlier. Finally, using the maximum 
principle and a continuity argument, it may be shown that there is a critical value of A, 
say A c > 0, such that there is no solution when A < A c and solution exists when A > A c . 
See [CCDiLtl] for details. Later, it was shown in [HD] that there is solution at A = A c as well. 
More results on existence and asymptotic behavior of doubly periodic solutions for (12.3p can 
be found in [2HE0j- Thus our understanding about the scalar case (12.2p or (12.3p is fairly 
satisfactory. 
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When G = SU (3) so that n = 2, things already become rather complicated because now 
one needs to resolve two coupled quadratic constraints. This more complicated problem was 
studied by Nolasco and Tarantello [51] who refined and improved the inequality-constrained 
minimization method developed in [10] for the n — 1 situation and showed that solutions 
in this n — 2 situation exist as well when A is sufficiently large. Note that, since in this 
case we are treating a system of equations, for which the maximum principle cannot be used 
and thus a continuity argument as that in [10] is not available, an existence result under the 
condition that A is sufficient large may be the best one can hope for. See PI for some new 
development regarding more generalized 2x2 systems arising in the Chern-Simons theory. 

The contribution of the present article is a successful settlement of the situation when the 
gauge group G is any compact group and in particular the Cartan matrix K of the equations 
is that of an arbitrary simple Lie algebra, of rank n. The system now consists of n nonlinear 
equations and results in n quadratic constraints, which cannot be resolved explicitly as 
hr [TUI 151] . In order to unveil the constraints difficulties we encounter and the varied levels 
of effeciency of the implicit constraint-resolution methods we use, we shall present a special 
method that works for the S77(4) situation, which has been extended to tackle the SU(N) 
situation in [34], and then a general method that works for all situations. The special 
method may be described as a “squeeze-to-the-middle” implicit-iteration strategy whose 
validity depends on the structure of the Cartan matrix of SU(N). The general method, on 
the other hand, uses a degree-theorem argument, which does not depend on the detailed 
specific numeric structures of the Cartan matrix. Rather, we shall see that, for a simple 
Lie algebra (say), things may be worked out miraculously to ensure an acquirement of all 
the needed apriori estimates so that the multiple quadratic constraints allow an implicit 
resolution in any situation under consideration. 

3 Chern—Simons—Higgs equations and existence theo¬ 
rem 

Our purpose is to carry out a complete resolution of the existence of doubly periodic solutions 
to (12.ip with very general Cartan matrix K. In order to treat the system in a unified 
framework, we need some suitable assumption on the matrix K. For a semi-simple Lie 
algebra, we know that the associated Cartan matrix has the property: the diagonal entries 
Ka assume the same positive integer 2, all off-diagonal entries Kij (i j- j) can only assume 
the non-positive integers —3, —2, —1, 0, and Kji = 0 if K t] = 0. This motivates us to consider 
(12. 1 p with a general matrix K, which satisfies 

K T = PS, (3.1) 

where P is a diagonal matrix with 

P = diagjPi,..., P n }, Pi> 0, % = 1,..., n, (3.2) 
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S' is a positive definite matrix of the form 

^ a ll ~ a 12 


S = 


—an —a{2 


Q-ii 



-OLi 


\ &nl &n2 . &nn— 1 ^nn ) 

OLn > 0, i — 1 ,..., n, a tj = > 0, i ^ j = 1 ,..., n, 

and all the entries of S^ 1 are positive. 


(3.3) 


(3.4) 

(3.5) 


In fact the assumptions on the matrix on K are broad enough to cover all simple Lie 
algebras realized as A n , B n , C n , D n , E 6 , E 7 , E 8 , F 4 , G 2 [MiEEQj thanks to the work of Lusztig 
and Tits [45]. 

As a result of (I3.5|) . the entries of ( K T )~ 1 = S^P^ 1 are all positive. In particular, we 
have 

n 

R i = y^((AT r )~ 1 )ij >0, i = 1,..., n. (3.6) 

3 =1 

Here is our main existence theorem for (12.11) . 


Theorem 3.1 Consider the non-Abelian Chern-Simons system (12.ip over a doubly periodic 
domain 12 with the matrix K satisfying (13.If) (13.51) . Let pa,...,p t N, L (i — 1,... ,n) be any 
given points on 12, which need not to be distinct. Then there hold the following conclusions. 

(i) (Necessary condition) If the system (12.ip admits a solution, then 


A > Aq = 




1B EE pP(k 

I67T i= 1 j= 1 


| 12 | 


E E \K 

i=l j=l 


(3.7) 


(ii) (Sufficient condition) There exists a large constant Ai > Ao such that when A > Ai 
the system (12.ip admits a solution (u),... ,«))). 

(iii) (Asymptotic behavior) The solution (u),... ,u)f) of (12.ip obtained in (ii) satisfies 


(■ e Ui — Ri) 2 dx —> 0 as A —> oo, i — 1,..., n, 


(3.8) 


Jn 

where R, (i — 1,... ,n) are defined by (13.6p . 

(iv) (Quantized integrals) The solution (u),... ,u(() of (12.11) obtained in (ii) possesses 
the following quantized integrals 



dx 


j =1 k =1 


47 rNi 
A 


i = 1,..., n. 


(3.9) 


In the subsequent sections we prove the theorem. 
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4 Necessary condition and variational formulation 


In this section we shall find a necessary condition for the existence of solutions of (12.1 ft and 
present a variational formulation. 

For convenience, we first use the translation 


UiUi + In Ri, i = l,...,n, (4.1) 

to recast the system (12.11) into a normalized form: 

( n n n \ Ni 

Y, K jk K^e Uk - Kije Uj + 4vr ^ 5 Pij (x), i = l,...,n (4.2) 

3= 1 k= 1 j= 1 / i=l 

whose vector version reads 


Au = \K\JK(XJ - 1) + 4vrs, 


(4.3) 

where the notation 

K = K T R = PSR, R = diag{i?i,..., R n ), 


(4.4) 

u = (m,... ,u n ) T , U = diag{e ui ,... ,e Un ), U = (e ui ,.. 

■ ■ i e“ n ) r , 

(4.5) 

/ M N n \ T 

1 — (1, • • • , 1) , S — I Spu, ■ ■ ■ > &Pns j 1 

\s=l 8 =1 / 


(4.6) 


will be observed throughout this work. Note that, since the matrix S is positive definite, so 
are the matrices 


A = p- 1 S~ 1 P- 1 and Q = RSR. (4.7) 

To solve (14.2ft or (14.31) over a doubly periodic domain, we need to introduce some back¬ 
ground functions to remove the Dirac source terms. Let u® be the solution of the following 
problem [3] 

= 47r^(5 Pis - J^u° i dx = 0, 

and Ui = u® + Vi, i = 1,... ,n. In the sequel we will use the n -vector notation 

V =(«!,..., u n ) T , N ={N 1 ,...,N n ) T , 0 = (0,..., 0 ) T . (4.8) 


Thus the system (14.2|) or (14.3j) becomes 


= A ( E E K jk K ije ^e<^ - ^ KijS 

<.3=1 k=1 3=1 


+Vj 


+ 


4tt Ni 

W 


i — 1,... ,n, (4.9) 















or 


Av = \KXJK(U - 1) + 


4vrN 

~W 


where 


U = diag{e u ° +vi ,..., e u - +Vn }, U = ( 


= (e u “ +v \...,e 


U°+ y n\T 


(4.10) 


(4.11) 


We now unveil a necessary condition for the existence of solutions of (12.11) . To this end, 
we rewrite (I4.10j) . after multiplying both sides of (14.101) by A, equivalently as 


AAv = AUQ(U-1) + —, 


where A and Q are defined in (14.7p . and 

b = (b u ..., b n ) T = 4vrAN = 47rP _1 5 _1 P _1 N. 
Noting (13.5j) . we obtain 

bi > 0, i — 1,..., n, 


(4.12) 

(4.13) 

(4.14) 


which will be used in the sequel. 

Let v be a solution of (14. lOj) . which is of course a solution of (14. 12 H . Then integrating 
(I4.12j) over O, we obtain the natural constraint 


which implies 


We may rewrite (I4.16P as 


/ UQ(U - l)dx + - = 0, 


f l T b 

/ U T Q(U - l)dx + —- = 0. 
'n * 


(4.15) 


(4.16) 


^ u -i)A( u -i| d * 


101 , 47Tl r p- 1 (A' r )- 1 N 

-—-l ?" 1 if _1 1- e --- 

4 v ' A 


where we have used the fact 


(iW)" 1 ! = = R1 


, (4.17) 


(4.18) 


and (14.13j) . 

Since the matrix Q is positive definite, (14. 1 7j) gives a necessary condition for the existence 
of solutions of (14.101) : 


101 lx , 4 7 rlLP- 1 (iW)- 1 N 

(/( )“ 1-4- L -> 0, 
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which establishes (i) in Theorem 13.11 

Now we show that the system (14.10)1 admits a variational formulation. To see this we 
consider the system (14.10)1 in its equivalent formulation, (14.12p . 

Now since the matrices A and Q defined in (1 4. 7jl are symmetric, we see that the equations 
(14.12)1 are the Euler-Lagrange equations of the functional 

2 

/(v) = [ diV T Adivdx + ^ [ (V - l) r Q(U - l)dx+ [ ^-^-dx. (4.19) 

^ (=] Jn * Jn Jn l“l 

Here and in what follows we use the notation (14. 7p . (14.8p . (14. lip and (14. 13p without explicit 
reference. 

We observe that the functional (14.19p is not bounded from below. So we cannot conduct a 
direct minimization. To deal with this problem, we will find a critical point of the functional 
/ by using a constrained minimization approach developed in pm, later refined by EH. 
Recently such a treatment was extended by [32] to solve the system associated with some 
general 2x2 Cartan matrices. To carry out this constrained minimization, the main difficulty 
is how to resolve the constraints, which will be the focus of the next three sections. 


5 The constraints 


In this section we identify a family of integral constraints under which our variational func¬ 
tional will be minimized. 

We start by decomposing the Sobolev space W 1,2 (Q) into W 1,2 (p) = M©H /1,2 (H), where 


W 1 ’ 2 (0) = w G W 1 ’ 2 (0) 


wdx = 0 


is a closed subspace of W 1,2 (h2). Then, for any Vi G VE 1,2 (J2), we have Vi = Ci + w^Ci G 
M, Wi G H /1,2 (h2),i = 1,... ,n. To save notation, in the sequel, we also interchangeably use 
W 1 ' 2 (H), W 1,2 (Q) to denote the spaces of both scalar and vector-valued functions. Hence, 
if v = w + c G W 1,2 (h2), with w = (wi,... ,w n ) T G H /1,2 (H) and c = (ci,..., c n ) T G 
satisfies the constraint (14 .15 p . we obtain 


diag{e ci ,... ,e Cn }Q 


( e Cl \ 


\e Cn ) 


- P 1 Rdiag{ai,... ,a n } 


/ e Cl \ 


\e Cn J 


A 


= 0 , 


(5.1) 
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where 


Q = Q(w) = RSR, 

( Quail — Ql2 a 12 


5 = 


— OLildi\ —Oii2di2 




^ln^ln 




\ 


(5.2) 


(5.3) 


\ C^nl^nl ^Ln2^n2 * * * ^nn^nn / 

and we adapt the notation 


di = aAwi) = / e“ i+w 


‘dx, 


dij = dij(wi,Wj)= / e“ i+ ^' +u,i+ ^da;, i,j = l,...,n. 

Jn 


Since the matrix S' is positive definite, from (15.21) and (j5.3[) we see that 

Q is positive definite. 

Now the system (15.ip can be rewritten in its component form: 


e ZCi Rfotada — e Ci ( ^ v Cj I!, I{j(\jj(ijj j + -A — 0, i — 1,..., n. 


Pi 




X 


(5.4) 

(5.5) 

(5.6) 

(5.7) 


For any w G kF 1,2 (fl), we see that the equations (15. 7 j) with respect to c are solvable only 


if 


Ridi 

IV 


E c D r> 1 ^ 4 R^biOliidn . . N 

e ] RiRjOiijdij \ > ---, i = l,...,n. (5.8) 




A 


In order to ensure (15.8ft . it is sufficient to take the following inequality-type constraints 


> 


Ary ■ ■ P^h 


A 


, i = 1,... ,n. 


(5.9) 


Define the admissible set 

= jw|w G W 1 ' 2 ^) such that (I5.9[) is satisfied j. (5.10) 

Therefore, for any w G srf, we can obtain a solution of (15. 7j) by solving the system 

1 I ( Ridi 


e Ci = 


2 Rj^ CXiidii | \ Pi 


Y^R'-rRR 


■j CX-ijCLij 




A 


Ridi \ -v „ „ \ 4 biRfaadi 

—p ^ / ; e 3 RjRjOijjdij j — 

* / 


= /?.(e Cl , • • •, e Cn ), i = l,...,n. (5.11) 

In the next two sections we aim at resolving the constraints (|5.7[) by solving (15. lip . 
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6 Resolving the SU( 4) constraints 

In this section we present a direct/concrete method for resolving the constraints (15.71) when K 
is the Cartan matrix of 577(4). Since in this case the coupling between the equations enjoys 
some special properties, we will see that the constraints allow a “squeeze-to-the-middle” 
solution process to be effectively carried out, which is of independent interest. 

For 517(4), the associated Cartan matrix K = ( Kij) is given by 



( 6 . 1 ) 


Obviously, now K satisfies all the requirement in Theorem 13.11 with P — I. Note that, in 
this case, 



( 6 . 2 ) 



(6.3) 



(6.4) 


Then the constraints (15. 7 j) are 



(6.5) 


( 6 . 6 ) 


(6.7) 


Furthermore, for any w G fF 1,2 (r2), to ensure the solvability of the equations (16.51) (16.71) 
with respect to c, the required inequality-type constraints (15. 9p take the form 



( 6 . 8 ) 


The admissible set srf reads 


srf = |w|w G W 1,2 {yi) such that (16.81) is satisfied j. 


(6.9) 
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Hence, for any w G &/, to get a solution of (j6.5j) - (16.7[) . it suffices to solve 


e Cl = 


e C2 = 


e C3 = 


«1 _|_ 2ai2.gC2 _|_ (— -\- 2ai2 c r °) 2 _ 8^1 a 22 

3 3 y V 3 3 / 9A 


2 a 


= A(e C2 ), 


li 


02 _|_ 3012 c Ci _|_ 3o23gC3 _|_ . / ( 22 - -J- 3a 12 QC1 _|_ 3023 C ct \ 1 _ 

48 8 V'48 8 / 2A 


2 a 


22 


03 2023 „C9 _i_ / 1 03 2 o 2 3 _ 863033 

3 1 " 3 e r VU r 3 e l 9A 


2a 


33 


= / 3 (e C2 )- 


( 6 . 10 ) 

= / 2 (e Cl ,e C3 ),( 6.11) 
(6.12) 


To solve these equations, we simply need to find a positive zero of the function 

F(t) = t - /a(/i(t),/ 3 (i)), i e [0, 00 ). 


(6.13) 


We have 

Proposition 6.1 For any w G as/, there exists a unique positive solution to for the equation 

F(t) = t-f 2 (h(t),Mt))=0. 

In view of this proposition we see that, for any w G the system fj6.10p ~ fj6.12p with 
respect to c admits a unique solution. However, as shown in the next section for the general 
case, there is no guarantee for the uniqueness of a solution to the general system (15.1111 . 
Proof. For any w G we see from (16.1011 - (16.1211 that 


/i(t) > 0, Vt >0, i — 1,3, / 2 (t,s) > 0, Vt, s > 0, 


which imply 

A direct computation gives 


Then, using the above limits and the Holder inequality, we have 


lim 


m 


t ~>00 t 


1 — lim 

t—> 00 




= 1 - 


°12 


l 23 


4aii a 22 4a 22 a 33 

1 1 

1 - 7 - 7 > °- 
4 4 


(6.14) 


no) = -/ 2 (/i(o),/ 3 (o))<o. 


lim /lW = 2 “ 12 , 
t—>c> 0 t 3&w 

(6.15) 

^ f 2 (t,t) 3(aq 2 + a 23 ) 

t^-oo t 8(222 

(6.16) 

Um / 3 « = 2o 23 
t 3a 33 

(6.17) 
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(6.18) 


































Then we see that 


(6.19) 


lim Fit) = oo. 

t—>oo 

Since we have F(0) < 0, then the function F(-) admits at least one zero to £ (0, oo). 
Now we prove that the zero of F(-) is also unique. In fact we easily check that 


d/i(0 



d t 



df 2 (t,s ) 



dt 

^(f + ^t + ^s) 2 - 

b2CL22 

2A 

df2(t,s) 

^Tf2(t,s) 


ds 


b2CL22 

2A 

dfsit) 

*Tf3{t) 



(£3 I 2 a 23 +\ 2 _ 8^3033 

V 3 -r 3 V 9A 


( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 

(6.23) 


which are all positive. Thus, the functions fi(t), (i = 1, 3) are strictly increasing for all t > 0. 
Then from (j6.20p (j6.23j) . (I6.10|l — (16.12j) and the constraints (16.8jl we obtain 


d m 

d t 


= 1 - 


> 






&2Q22 

2A 


X 


“ 12/1 (*) 


LUf + M _ 


2ai2 +\ 2 _ 8&ian 

9A 


: + 


a; 


'23 


Mt) 


v/(f+ ¥ 0 2 - 


8&3 a33 

9A J 


> 1 




a 12fl(t) + ^23 f3^) 


^2^22 

2A 


Ut + MM + ttAW) 

^ _ |/2(/i(t), fsjt)) anfijt) + «23/3(^) 

^fl(t) + *f-f 3 (t) - 


= 1 - 


/ 2 (/i(t),/ 3 (t)) F(t) 


t > 0, 


so the uniqueness of the zero of F(t) over [0, oo) follows from the monotonicity of F(t)/t for 
t > 0 and the proof of the proposition is complete. 

Using Proposition 16.11 for any w £ stf , we see that the equations (16.51) 116.71) with respect 
to c admit a unique solution c(w) = (ci(w), c 2 (w), c 3 (w)) r determined by (l6.10jMI6.12D . 
such that v = (ui,U2,u 3 ) r defined by 


Vi = Wi + Cj(w), z = 1,2, 3, 


(6.24) 


satisfies the constraints (14.151) for the SU (4) case. 
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7 Solving the constraints in general 


In this section we carry out a new way to resolve the constraints (15.7ft by solving (15. lift 
for the general case via a topological-degree-theory argument. To this end, we consider the 
system 

F(t) = t — f(t) = 0, t = (t u . .. ,t n ) T G K", (7.1) 

where = (M+) n , f(t) = (/i(t),. .., f n ( t)) r . We have 

Proposition 7.1 For any w e , the system (17. Tj) admits a solution t e (0, oo) n . 

Proof. To conduct a degree-theory argument, we deform the system (17.ip as 

F(e,t) =t-f(e,t) = 0, t G M+, ee[0,l], (7.2) 

where 


f(e,t) = (/i(e,t),...,/ n (e, t)) T , 
/i(e,t) = 


2 Rj, OLH&H 


I ^ + v tjRtRjai , 0i . j 


A 


/ v tj-Lli-TljO'ij 


Pi 




AebiBjauaa 

A 


i — 1 ,..., n. 

Then, we see that, to solve the system (17.1ft . we need to find a solution of 

F(l, t) = 0, t e Rj. 

To facilitate our statement, we use the convention that we write 


(«i, ..., a n ) T < (<) (/3i, ..., f3 n ) T if ai <{<)Pi,i = 1, — , n, 

and we use the same convention for matrices. 

To proceed, we establish the following key a priori estimates. 

Lemma 7.1 For any w e stf and e G [0,1], every solution t of (17.2p satisfies 

0 < aJi < |0|, i — 

0 < U < 1, i — 1 ,..., n. 


By virtue of Lemma [7.11 we immediately obtain the following. 


(7.3) 


(7.4) 


(7.5) 


(7.6) 

(7.7) 
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Corollary 7.1 For any w G s4, every solution (ci,..., c n ) of (15. lip satisfies 

aie Ci < |f2|, i — 1,..., n, 
e Ci < 1, i — 1,..., n. 


(7.8) 

(7.9) 


Proof of Lemma 17.11 For any w G and e G [0,1], let t be a solution of (17.2ft . By 
(17.21) (17.41) we readily get the left-hand sides of (I7.6j) — (j7.7p . 

From (17.2p it is straightforward to see that 


pF + I] 


U = fi(e,t) < 
which can be rewritten as 




Rj O-H On 




R t auaati ^ ^ tj Rj Rj riy — 




Ri&i 




(7.10) 


(7.11) 


Noting the expression of Q, we rewrite (17. lip in a vector form 

Qt < P~ 1 R&, 

where a = (ai,..., a n ) T . By Holder’s inequality we obtain 

1 1 

<% < OiiQjj, ", £ -W i,j = l,...,n. 


(7.12) 


(7.13) 


Therefore, noting that tj > 0, i — 1,..., n, the expression of Q, and (17.13p . we arrive at 


(7.14) 


diag ja^,..., a^ n | Qdiag ja^,..., a£„j1 < Qt. 

On the other hand, by (13.51) and the expression of Q in (14.7p we see that 

all the entries of <2 _1 are positive. 

Therefore, by (I7.12p . (I7.14p . and (17.15j) . we get 


(7.15) 


t < 


diag { a n 2 ,..., a n n ) Q Miag { a u 2 ,..., a n f } P 1 R& 


= diag |a n 2 ,..., a n n j Q~ diag aia n 2 ,..., a n a n n \ P~ R 1. 
From (17.13p . (I7.15p . and (I7.16P we infer that 
diag{ai,..., a„}t 


(7.16) 


< diag jaia n 2 ,..., a n a n f \ Q diag \ aia n 2 ,..., a n a n f \ P R1 

< {niQ-'p-'Ri 

= |n|i, 


(7.17) 
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where (14.18p is also used . Then the right-hand sides of (I7.6[) follow from (17.171) . By Jensen’s 
inequality, we have a* > |12|,i = 1 ,,n, which together with (17.171) imply the right-hand 
sides of (17.7p . The proof of Lemma mi is complete. 

By the definition of F(e,t), for any e G [0,1], it is easy to see that F(e,t) is a smooth 
function from R” into R”. 

Let us define 

fU(0,r o r, (7.18) 

where ro > 1 is a constant. 

Then by (17.71) we see that, for every w G stf and £ G [0,1], F(e,t) has no zero on the 
boundary of 12. Consequently, the Brouwer degree 

deg(F(e, t), 0, 0) (7.19) 

is well dehned. 

To prove Proposition 17.11 it is suffhcient to show that 

deg(F(l,t),fi,0)^0. (7.20) 


Since F(e, t) is smooth, by homotopy invariance [T8], we have 

deg(F(l, t), O, 0) = deg(F(0, t), 12, 0). (7.21) 

Now we only need to calculate the degree deg(F(0, t), 12, 0). 

Note for any w G srf the system 

F(0, t) = 0 (7.22) 


is reduced into 


U 


~pT + X/ tjRiRjOiijCLij 


= 0, i — 1,..., n. 


(7.23) 


Ri 

By the definition of Q we rewrite the system (I7.23P equivalently in a vector form 

Qt = P~ l Ra. (7.24) 

In view of (15.6p the matrix Q is of course invertible. Then we see that the system (I7.24p . 
i.e., (I7.22p . has a unique solution 


t = Q-'p-iRa, 


which, belonging to 12, is not a boundary point of 12 by (17.7j) . 

Noting (15. 6p . the determinant of Q is positive, which implies that the Jacobian of 
F(0,t) is positive everywhere. Therefore, by the definition of the Brouwer degree, we have 
deg(F(0, t), 12, 0) = 1, which implies deg(F(l, t), 12, 0) = 1. Then the proof of Proposition 
17. 1 1 is complete. 

By Proposition 17.11 we see that for any w G the constraint equations (15.71) admit 
a solution c(w) = (ci(w),... ,c n (w)) T determined by (15. lip , such that v = w + c(w) = 
(wi + Ci(w),..., w n + c n (w)) T satisfies the constraints (I4.15p . 
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8 Constrained minimization 


In this section we solve the equation (14.9ft by finding a critical point of the functional / via 
a constrained minimization procedure. To do this, we consider the constrained functional 


J(w) = /(w + c(w)), w G sA , 


( 8 . 1 ) 


where c(w) is the solution of the constraint equations (15.7ft determined by Proposition 17. 11 
Since, for every w G srf, v = w + c(w) satisfies the constraints (14.151) . we have 


(U - 1) T Q(U - l)cto = / 1 T Q(1 — U)da: — 


l r b 


I l T P” 1 i?(l — U)dx 

Jn 


l r b 


Then the functional J can be expressed as 

2 


J(w) = - V' f diW r Adiwdx+—l T P ' R f (1 — U)dx + b T c 
2 . =1 Jn 2 J n 


rb 


\tS a ^ rAB ^ ix+ \t}jA 


<Ci Ruy-\-w. 

e e ’■ 


; )dx 


J2 hiCi ~ 9^ bi ’ 


1=1 


i =1 


( 8 . 2 ) 


(8,3) 


where and in the following we use the notation (14.71) . (14.13[) . 

We easily see that the functional J is Frechet differentiable in . We aim to show that 
the functional J admits a minimizer, say w, in the interior of srf. Then v = w + c(w) is a 
critical point of the functional /. 

To proceed further, we need the following inequalities. 


Lemma 8.1 For any w e and s G (0,1), there hold the inequalities 


U^+Wi 


dx < 



i — 1 ,..., n. 


Such type of inequalities were first established in [49] . 

For our purposes, we will also need the Moser-Trudinger inequality [26] 


(8.4) 


e w dx < C exp 



Vm e 


(8.5) 


where C > 0 is a constant. 

Let a 0 and (3 0 be the smallest eigenvalues of A and Q defined by (14.71) . respectively. 
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Lemma 8.2 For every w e srf, the functional J satisfies 

n 

J( w )> y^llV^II^-COnA+l), (8.6) 

i=l 

where C > 0 is a constant independent of A. 


Proof. By (|4.7[) and the definition of J we see that 


</( w ) > + 


i =1 


i=l 


Noting (15. lip we obtain 

e Ci > 

which together with (j5.9|) yield 

c . ^ 2 PA 
e 1 > 


2 1 ( Hj o 77 a {{ 


, i = 1 ,... ,n, 


2 PA 


XRiOi XRi f Q e u i +Wi dx 


, i = 1,... ,ra. 


Then we have 


Ci > In 


2 PA 
Ri 


In A — In / e Ui+Wi dx, i — 1,..., n. 


Using Lemma [8. II and the Moser-Trudinger inequality (18. 5p we arrive at 

1 
s 


In [ e u < +Wi dx < -—- 


{In A - In (4Pj l b i a a ) } + - In / e su ° +SWi dx 


< 


-|-||V^||| + — {inA — In (4 Pfaaa)} 
167T S 


In C 


S f2 

Combining (18. 7ji . (18.1 Oh . and (18. lip , we have 


+ maxtq, i — 1 ,..., n. 


(8.7) 


( 8 . 8 ) 


(8.9) 


( 8 . 10 ) 


( 8 . 11 ) 


X w) > I 21 


S max{6j}\ n n 

i<i< n \ ^ || V Wj ||2 - - ( ln A - ln ( 4P iVn) + InC) 

i =i s i =i 


2 16vr 

n 

y bi jin ( 2R i P i b i au ) + 


2=1 


max u, 
n 1 


( 8 . 12 ) 


Therefore, taking s suitably small in (18.121) we obtain the lemma. 

Using Lemma 18.21 we see that the functional J is bounded from below and coercive in 
sF. Then noting that J is weakly lower semicontinuous in &/, we conclude that J admits a 
minimizer in sF. In the following we establish some estimates as in [T01I5T] to show that this 
minimizer belongs to the interior of srf when A is sufficiently large. 
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Lemma 8.3 There exists a constant C > 0 independent of A such that 


|S1|A . / R, 

mi J w) > - mm < — 

wG&z/ 2 l<i<n I P; 


-V l — C(1 + In A + VX). 


(8.13) 


Proof. By the definition of &/, we see that at least one of the following equalities 

(8.14) 


of 4 aaPfhi . 

— =-t-, i = n 

da A 


must hold on the boundary dsrf. 

If the case i — 1 in (18. 14[) holds, by (17.171) we obtain 


TD TD 

aie ci < + ^2-^-(Q- 1 ) lj a 1 a j 


1 _ 1 
2. 2 


< 


Pi 

R\ 
~Pl 


Pj 
3= 2 3 


a u a jj 


n TD 

, 1 X-> Ibn 


( Q 1 )n«?a n 1 +|n |2 1 )ijOia 11 ' 

i=2 1] 


= (Q ) 


_in 4Pii2i6ian 


U D 

2 PiE^(Q" 1 ; 


ii" 


A 


+ 


3 =2 




7a 




If other cases happen, similar estimates as (18. 15j) can be established. 
Using (17. 8p and (18. 15[) leads to 


(8.15) 


A 

2 


V ^ [ (1 - e Ci e u * +Wi )dx 

i=i ^ 



(8.16) 


By (18.16[) and estimating q (f = 1,..., n) as that in Lemma [8.21 we get the lemma. 

Now we need to choose some suitable functions in the interior of stf to estimate the value 
of the functional J in srf. 

Recall that in [60j Tarantello proved that, for /j, sufficiently large, the problems 

Av = M e“?+”(e""+” - 1) + i = (8.17) 

admit solutions vf{i = 1,..., n), such that vf + vf <0 in Q, cf — Ur j n vfdx 0 and 
wf — vf — df —y —pointwise as/i—»oo,i = l,...,n. Then it follows that 

lim [ e u °i +w idx = |0|, lim [ e u * +u< i +w * +w j dx = |f2|, i,j = l,...,n. (8.18) 

J n M->°° J n 

Consequently, in view of the expression of Q , we have 


lim Q( w M ) = |f2|Q. (8.19) 

fl—> OO 
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( 8 . 20 ) 


Noting (15.61) . we know that Q( w M ) is invertible, which together with (j8. 19j) yield 


lirn Q l (yr tl ) = —Q '. 

11 11 

Therefore, by the definition of srf, (18. 18ft and (18.201) . we conclude that, for a fixed Aq > 0 
sufficiently large and any e e (0,1), there exists a fj, e 1, such that 


w Me = (u>i £ ,..., iu£ e ) T G int^, 


for every A > Ah and 


«yK £ ,™f) < (l + e)|fi| < 2|fi|, i,j = l,...,n, 
( 1_£ )^-l ^ „ ( 1+£ )^-l ^ 2 r»-l 


|n| 


-g- 1 < g-^w^) < 


VI ‘ e < N ' Q 


( 8 . 21 ) 

( 8 . 22 ) 

(8.23) 


At this point we can establish the following comparison result . 
Lemma 8.4 For w Me given by (18.211) . we have 


J(w^)- inf J(w)<-1, 

w E.O'W 


(8.24) 


when A is large enough. 


Proof. Since w Me e int^, by (15. lip and the Jensen inequality, we arrive at the 
estimates 


^ + £ e c ^ s) RiR 


3 C<( W^) _ 


3+i 


j Q-ij &ij 


OiiiCLii 


t 


X 


\ 


1 + 


1 


1 - 


Rj OtiiQ’ii 


A ( ^pr- + y o'' R; IIj(\ijai, 

jA* 


Pi 


> 


VV'.< " Il.IljOjjUij 
j¥=i _ 

Ri 


2b, 


A -^r 1 + XI e Cj RiRjOtijOij 

1 j# 


Pi 


> 


+ E e^^RiRjOijOij 
j¥=i 


2 PA 

A|n|i2, 


■, i = 1,..., n. 


Here and to the end of this section we understand that 

Q>i = a i{ w i S )i aij = a ij( w i e J w j S )i i,j = 1,... ,77,. 


(8.25) 


(8.26) 
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Then it follows from (18.25p and (j8.22[) that 




> 


Pi A|ft| 

Ri\fl\ AoiaPiRibi 
P 


A 


, i = l,...,ra,(8.27) 


which can be expressed in a vector form: 


Q{ w*)(e' 


ftV 0 ci( w ^)^ _ ^ e c n (w ^)\r 


> ^(P^Pl 


4PP 


A 


diag {an 5 • • • j ^rm } b. (8.28) 


Noting that all the entries of Q 1 are positive, using (I8.28j) and (18.231) . we have 


( e c i( wME ) gCnlw^nr 


> (Qig-^w^^p^pi 


AQ-\w^)PR 


A 

8 Q-'PR 


diag {an 5 • • • ? &nn} b 


> (1 - e)Q 1 P 1 R1 - - diag {an, ■ • •, a nn } b 


A|n| 


= (1 -e)l 


8 Q-'PR 


A|Q| 


diag {an •)••••> &nn} t). 


(8.29) 


Hence from (18.291) we get 

n Q 77. 

/ (l - e Ct(w "' £) e u > +w '^ dx < \n\e + j ^(Q-^P.P^a^, * = 1, ■ ■ •, n. (8.30) 


3 = 1 


Using (17.9j) and (18.3011 . we see that there exists a constant C £ depending only on e such 


that 


+c. 

3 = 1 1 


(8.31) 


Then it follows from Lemma 18.31 and (18.311) that 


•/(w“-)- mf J(w)<12h 

W z 


\ ^ Ri 

}^—£- mm 


o'=i 


Pi l<i<n 1 p 


Pi 


+ C (\/A + lnA + l) , (8.32) 


where C > 0 is a constant independent of A. Now taking £ suitably small and A sufficiently 
large in (18.321) . we get the lemma. 

Therefore by Lemma 18.21 and Lemma 18.41 we conclude that there exists a large Ai > 
max{Ao, A' 0 } such that, for all A > Ai, the functional J has a minimizer in the interior of 
say 

w A G kite/. (8.33) 

ft is straightforward to check that 

v A = w A + c(w A ) (8.34) 

is a critical point of / and accordingly a solution of the system (14.101) . Hence the second 
conclusion of Theorem 13.11 follows. 
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9 Asymptotic behavior and quantized integrals 


In this section we prove the last two conclusions of Theorem 13.11 We first establish the 
asymptotic behavior of the solution obtained above as A —* oo. 

Lemma 9.1 Let v A be given by (18.34)1 . Then, 


lim 

A—>-oo 


( e u i+ v t 


l) 2 dx = 0, i — 1,..., n. 


(9.1) 


Proof. For any e G (0,1), we conclude from (18.311) that there exist constants A e > 0 and 

\n\Xe ^ Ri 


C e > 0 such that 


J(w A ) = inf J(w) < 

Z L ' K 

i=i 1 


£G+o, 


for all A > A e . 

Since Q is positive definite and /3 0 being the smallest eigenvalue of Q, we get 
f(U- 1) T Q(U - l)cto > p 0 V f {e u ° +Vi - l) 2 da;. 

J f2 j = j J f2 

Estimating c t as that in Lemma [8721 and using (19.31) . we see that 

J(w A ) > — V [ (e“°+^ A - l) 2 da; - C7(ln A + 1) 

where C > 0 is a constant independent of A. 

Therefore, we infer from (19. 2 p and (19.4[) that 

lim sup ^ [ (e u * +v i — l) 2 dx < p, Ve G (0,1). 

\-+oo SI ~ Po tl ^ 


(9.2) 


(9.3) 


(9.4) 


(9.5) 


Since £ G (0,1) is arbitrary, the lemma follows immediately. 

Then using the translation (14.ip and Lemma 19.11 we get the third conclusion of Theorem 

EH 

Finally, we can establish the quantized integrals (13. 9 j) . In fact, for the obtained solution, 
integrating the equations (14.91) . we see that desired quantized integrals follow. 

The proof of Theorem 13. II is now complete. 
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